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Abstract 

A  design  of  fuel-cell  manifold  to  evenly  distribute  the  reactants  among  the  channels  within  bipolar  plates  is  presented.  Neglecting  the  friction 
term  from  the  momentum  equation,  one-dimensional  governing  equations  are  manipulated  to  develop  the  manifold  equation  to  achieve  the  required 
design.  Some  manifold  shapes  that  satisfy  the  manifold  equation,  are  suggested,  together  with  the  distribution  of  static  pressure  within  the  manifold 
itself. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

A  fuel-cell  stack  has  a  hydraulic  network  to  distribute  the 
reactants  among  the  component  cells.  The  interface  between  the 
external  source  of  reactants  and  the  cell  inlet  and  that  between  the 
cell  outlet  and  the  exhaust  plumbing  are  each  known  generically 
as  the  fuel-cell  manifold  or  header.  Two  typical  configurations 
of  the  fuel-cell  manifold,  e.g.,  Z-type  and  U-type,  that  provide 
different  flow  directions  are  shown  in  Fig.  1.  The  research  to 
date  has  focused  on  developing  the  channel  distribution  for  a 
given  manifold  design  of  rectangular  cross-section.  This  study 
reported  here  shows  how  to  change  the  manifold  shape  to  achieve 
an  even  distribution  of  reactants  among  the  channels. 


reactants  while  the  lower  manifold  collects  them  after  they  have 
passed  through  the  channels.  The  governing  equations  for  the 
manifold  flow,  including  the  area  change,  can  be  written  as  fol¬ 
lows: 


Upper  manifold 


d  (bh  i  V\ ) 

dxj 


N 

L 


AcVc 


d  (pblnVf) 
d.vi 


d(  P\  bh  i ) 
d.ri 


fpVfo  +  hi) 


Lower  manifold 


(1) 

(2) 


2.  Mathematical  model 

2.1.  Z-type 

Maharudrayya  et  al.  [1]  have  proposed  an  analytical  solution 
for  manifold  flow  with  a  uniform  cross-section.  In  the  work 
presented  here,  the  manifold  is  allowed  to  change  its  cross- 
section  so  as  to  achieve  even  distribution  among  the  channels. 
A  schematic  diagram  of  the  proposed  system  and  its  key  dimen¬ 
sions  are  given  in  Fig.  2.  The  upper  manifold  provides  the 


d(bh2V2 ) 
dx2 


N 


TAcVc 


d(pb/i2V|) 

dx2 


d(P2bh2) 

d.V2 


fpV2(b  +  h2) 


(3) 

(4) 


For  convenience,  the  inlet  of  the  upper  manifold  is  aligned  with 
the  first  channel  and  the  outlet  of  the  lower  manifold  with  the 
last  channel.  The  consumption  of  mass  by  the  electrochemical 
reaction  is  assumed  not  to  affect  the  channel  flow.  The  pressure 
gradient  and  flow  velocity  is  given  by: 
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Pi-P2  =  AP12  =  H-pVf, 


where  H  = 


4  fLc 

Dc 


(5) 
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Nomenclature 

A  cross-sectional  area 

b  manifold  width 

D  hydraulic  diameter 

/  friction  factor 

h  manifold  height 

H  loss  coefficient 

L  manifold  length 

M  flow  parameter  defined  in  Eqs.  (8)  and  (9) 
N  number  of  channels 

P  static  pressure 

q  flow  rate  per  unit  width  of  manifold 

Re  Reynolds  number 

V  velocity  magnitude 

W  channel  width 

x  co-ordinate  along  the  manifold 

Greek  letters 

a  manifold  height  ratio 

(J>  manifold  height  ratio 

y  manifold  height  ratio 

A  difference  operator 

0  non-dimensionalized  channel  pressure  drop 

\x  viscosity 

77  non-dimensionalized  dynamic  pressure 

p  density 

P2  channel  pressure  drop 

Subscripts 

1  upper  manifold 

2  lower  manifold 

c  channel 


Fig.  2.  Schematic  diagram  of  manifold  shape  with  key  dimensions:  (a)  upper 
manifold;  (b)  lower  manifold. 


the  friction  term  can  be  compared  in  terms  of  order  of  magnitude 
by: 

friction  term  fpV2(b  +  h)  /  2 /L\ 

inertia  term  A(pbh  V1)/  Ax  \  ) 

Since  the  friction  factor  for  turbulent  tube  flow  is  below  0.01 ,  the 
friction  term  becomes  non-negligible  if  the  hydraulic  diameter 
of  the  manifold  is  reduced  to  below  2%  of  its  length.  In  this  study, 
the  manifold  is  assumed  to  have  a  cross-section  of  moderate  size 
so  as  to  cause  negligible  friction.  With  the  introduction  of  this 
assumption,  the  momentum  equation  can  be  simplified  to: 

Upper  manifold: 

q\ 

p—  +  P]  h  i  =  const  =  Mj  (8) 

h  i 


Vc  =  AP12 


Dr. 


2  pLc(Ref)c 


where  (Ref)c 


=  13.84+  10.38  exp 


(6) 


The  pressure  gradient  is  established  within  each  manifold  due  to 
the  inertia  change  and  the  friction  with  the  wall.  The  inertia  and 


Lower  manifold: 

q i 

p—  +  7*2^2  =  const  =  M2  (9) 

h  2 

To  facilitate  the  manipulation,  the  main  variable  was  changed 
from  the  velocity  to  the  flow  rate  per  unit  width,  q.  Eqs.  (8)  and 
(9)  are  used  to  calculate  the  driving  force  at  the  channel  or  the 


(a) 

Fig.  1.  Fuel-cell  manifold:  (a)  Z-type;  (b)  U-type. 
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pressure  difference  between  the  upper  and  lower  manifolds: 


AP12  = 


Mi 

7*7 


(10) 


The  continuity  equations  at  both  manifolds  give  rise  to  the  fol¬ 
lowing  expression  that  greatly  simplifies  the  calculation: 


Since  the  momentum  equation  does  not  change,  Eq.  (10)  can 
still  be  used  to  represent  the  pressure  difference  between  two 
manifolds  of  U-type  configuration.  The  continuity  equation  pro¬ 
duces  a  different  relationship  between  the  q\  and  qj  flow  rates, 
i.e.: 


q\=qi  =  q 


(19) 


qi  q2  —  const  =  q 


(11) 


The  flow  rate  q  has  the  value  of  q\  at  the  inlet  of  the  upper 
manifold  or  at  .v  =  0  since  qi  is  zero  there.  Using  Eqs.  (6)  and 
(10),  Eq.  ( 1 )  can  be  written  in  terms  of  a  single  unknown  variable 
q\  if  Mi  and  Mi  are  known,  i.e.: 


d£i  _  N  DgAc 
dx  bL  2 pJLc(Ref)c 

Ml  _  Ml  _  ( q\_  (q  -  gi)2^ 

hi  h 2  P  h\  J 


(12) 


It  is  worthwhile  nothing  that  q\  and  P\  are  not  given  simulta¬ 
neously  as  input  for  the  simulation.  Hence,  instead  of  using  Eq. 
(8)  directly  to  determine  Mi ,  it  has  been  decided  to  employ  the 
pressure  difference  at  the  end  of  the  upper  manifold  (or  at  x  =  L) 
and  the  reference  pressure  at  the  lower  manifold,  i.e.: 


Ml  =  /u,*=L(Pamb  +  (AP12W)  (13) 

Since  the  pressure  difference  is  usually  available  after  the  cal¬ 
culation,  the  final  solution  is  obtained  after  the  iteration.  Mi  can 
be  calculated  from  qj  and  Pi  at  the  outlet  of  the  lower  manifold 
where  these  two  parameters  are  the  flow  rate  at  the  inlet  of  the 
upper  manifold  and  the  ambient  pressure,  respectively: 

q2 

M2  =  P~ - h  Pamb^2,x=L  (14) 

hi_,x=L 


2.2.  U-type  manifold 


A  similar  formulation  to  that  above  is  possible  for  the  U-type 
configuration,  it  is  necessary  only  to  make  a  slight  modification 
due  to  the  flow  direction  at  the  lower  manifold.  For  convenience, 
the  inlet  and  the  outlet  are  aligned  with  the  first  channel  from 
the  left  in  Fig.  1. 


Using  Eqs.  (6)  and  (10),  Eq.  (15)  can  be  written  in  terms  of 
single  unknown  variable  q  when  M\  and  Mi  are  known,  i.e.: 

d q_  N  D2cAc  [Mi  Mi  _2  /  1  lY 
dx  bL  2pLc(Ref)c  hi  Y  ^  \h\  h\) 

(20) 


Mi  and  Mi  are  calculated  in  different  way  from  Eqs.  (13)  and 
(14),  since  in  the  outlet  of  the  lower  manifold  the  static  pressure 
is  now  located  at  x  =  0.  Thus: 


Mi—p 


M2  =  p 


q 


h  i,*=o 

,2 


q 


hi,x=o 


T  h  I  ,x=o(  Pamb  "T  ( A  /J|  2  )a=o) 


T  Y.a=0  Pamb 


(21) 

(22) 


3.  Optimization  of  manifold  shape 

A  main  goal  in  the  design  of  a  manifold  is  to  distribute  the 
reactants  among  the  channels  as  even  as  possible  with  a  moderate 
size  of  manifold.  As  the  manifold  size  increases,  it  will  certainly 
improve  the  distribution  but  will  become  useless  from  a  prac¬ 
tical  point  of  view.  By  contrast,  if  its  size  becomes  too  small, 
the  pressure  drop  will  be  raised  and  the  system  will  require  an 
uneconomical  air  blower  for  the  distribution  of  reactants  among 
the  channels.  In  this  section,  the  manifold  shape  is  optimized  to 
achieve  the  best  compromise  in  size. 

If  the  manifold  is  designed  to  distribute  the  same  amount  of 
reactant  among  each  channel,  then  the  channel  velocity  and  the 
flow  rate  through  the  manifold  can  be  written  as: 


Using  Eq.  (6),  the  pressure  drop  in  the  channel  is  determined 
without  iteration  as: 


Upper  manifold: 


d \qi 
dxi 


N 

bL 


ACVC 


q  7 

p—  +  Pi  hi  =  const  =  Mi 
hi 


(15) 

(16) 


D\  NAC 
3.1.  Z-type  configuration 

Substituting  Eq.  (23)  for  q \  in  Eq.  (12)  yields: 


(24) 


Lower  manifold: 


dc/2 

dxx 


N 

bL 


AcVc 


qi 

p-T1-  +  P2I12  =  const  =  M2 
h  2 


=  L2,  where  x  —  — 

L 

(25) 

The  manifold  length  L  was  used  to  normalize  the  stream-wise 
co-ordinate.  Rearrangement  of  the  terms  after  substitution  of  M\ 


(17) 


Mi 

~hl 


M2 

Il2 


-\2 


pq 


(i  -  x) 

,  hi 
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Fig.  3.  Profile  of  lower  manifold  for  h\  (x)  =  const  and  normalized  static  pressure  Fig.  4.  Profile  of  upper  manifold  for  h2(x)  =  const  and  normalized  static  pressure 

at  upper  manifold  (Z-type).  at  upper  manifold  (Z-type). 


and  M2  in  Eqs.  (13)  and  (14)  using  Eq.  (24)  gives: 


hl,x=l 

hi 


(^2  +  ^amb) 


hl,x=\ 

hi 


Pamb 


-pq 


2 


(  0  ~  x)2 

V  h\ 


(26) 


The  various  combinations  of  h\  and  hi  that  satisfy  Eq.  (26)  will 
render  an  even  distribution  among  the  channels.  From  the  defi¬ 
nition  of  M i ,  the  static  pressure  at  the  upper  manifold  becomes: 

P\  =  ( p amb  +  Q)  -  pCL(\  _  i}2  (27) 

hi  hi 


3.1.1.  hj(x)  =  const 

To  find  the  specific  design  that  satisfies  Eq.  (26),  the  special 

case  is  assumed  in  which  the  height  of  the  upper  manifold  is 

maintained  at  a  constant  value.  Then: 

2 

pq 

hi,x=i  =  ahi,  — - -x, - =  77, 

"amb  «  2, 1=1 


£2 

Farnb 


=  0, 


hi 

hi.x=l 


=  hi, 


(1  -  cr77(l  -  x)2)h\  -  (77  +  1)/j2  +  77x2  =  0  (28) 


During  the  derivation  of  Eq.  (28),  three  scaling  variables  are 
suggested:  (i)  a  is  the  ratio  between  the  heights  of  both  manifolds 
based  on  the  height  at  the  exit  of  the  lower  manifold;  (ii)  77  is 
the  ratio  of  dynamic  pressure  to  static  pressure  at  the  outlet  of 
the  lower  manifold;  (iii)  ©  is  the  pressure  drop  in  the  channel 
normalized  with  respect  to  ambient  pressure. 

The  static  pressure  at  the  upper  manifold,  normalized  to  ambi¬ 
ent  pressure,  is: 

7>i  =  1  +  0  -  a2  77(1  -  x)2  (29) 


The  height  profile  and  the  static  pressure  for  various  a  values  are 
plotted  in  Fig.  3  with  hi  fixed.  The  static  pressure  increases  along 
the  manifold  co-ordinate  x  to  maintain  force  balance  with  the 


reduced  inertia.  The  size  of  the  arrow  drawn  downward  between 
two  curves,  i.e.,  the  channel  pressure  drop  (72),  is  shown  to  be 
constant  along  the  manifold  co-ordinate  x.  To  achieve  even  dis¬ 
tribution  for  the  upper  manifold  with  a  uniform  cross-section, 
the  lower  manifold  was  modified  to  have  its  cross-sectional  area 
decreasing  along  the  co-ordinate.  If  a  decreases  or  the  size  of 
the  lower  manifold  is  reduced,  the  required  variation  in  shape 
becomes  large  in  order  to  provide  an  even  distribution.  If  the 
number  of  channels  or  their  area  is  changed,  the  scaling  vari¬ 
ables  can  be  used  to  rescale  the  manifold. 


3.1.2.  Ii2( x)  =  const 

In  this  case,  the  lower  manifold  has  a  uniform  area  while  the 
upper  manifold  is  allowed  to  change  its  cross-section  to  enable 
the  channel  distribution  to  be  even,  i.e.: 


h  i,x=i  = 


pq 


Q 

^amb 

(1  +  0 


=  0, 

77(1 


P'Mubhj 
hi 


=  77, 


h  l  ,x=i 
■  x2))h2 


=  h  i 


(1 


@)hi  +  /J277(l  —  x)2  =  0 


(30) 


In  this  relationship,  the  scaling  variables  are  defined  as:  (i)  j3 
is  the  ratio  between  the  heights  of  both  manifolds  based  on  the 
height  at  the  inlet  of  the  upper  manifold;  (ii)  77  and  0  have 
the  same  physical  meaning  as  in  Section  3.1.1,  but  with  a  dif¬ 
ferent  length  scale.  The  normalized  static  pressure  at  the  upper 
manifold  is: 

Pi  =  J-(l  +  0)  -  J=/i277(l  _  x)2  =  1  +  0  +  77(1  -  x2) 
hi  h\ 

(31) 


The  static  pressure  decreases  along  the  manifold  co-ordinate  x, 
as  shown  in  Fig.  4.  This  may  be  attributed  to  the  fact  that  the 
upper  manifold  is  designed  to  have  the  cross-section  increasing 
along  the  x  co-ordinate  and  the  exit  of  the  lower  manifold  is 
located  at  x  =  L.  As  will  be  shown  later,  the  shape  of  manifold 
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exerts  less  influence  on  the  distribution  of  static  pressure  in  a 
U-type  configuration  where  the  exit  is  located  in  the  same  co¬ 
ordinate  with  the  inlet  at  x  =  0.  The  shape  of  the  upper  manifold 
is  consistent  with  that  of  the  lower  manifold  in  Section  3.1.1 
because  of  the  symmetry. 


3.1.3.  hi(x)  =  h2(L  —  x) 

A  symmetrical  manifold  design  h\(x)  =  h2(L  —  x)  is  an  obvi¬ 
ous  consideration  given  the  location  of  the  inlet  and  the  outlet. 
Introducing  the  new  co-ordinate  x'  =  1  —  x  and  using  the  rela¬ 
tions  hi(x!)  —  /z 2 ( 1  —  x')  and  h 2(x')  =  h\(\  —  x'),  Eq.  (26)  can 
be  rewritten  after  dropping  the  prime  notation  as: 


h\  r 


l,x=l 


hi{x) 


^2,1—1 
h  i  (x) 


hi-x=\ 


Pamb 


X2 


-  pq 


d-x)2 


h2{x)  hi(x¥ 
Eq.  (26)  and  (32)  are  added  to  produce: 


=  Q  (32) 


h\,x=  t 

h  i  ,i=o 

h[ 


=  Y > 


=  hi 


pq~ 


=  n, 


Q 


T’amb^  T,i=0 

1  1 


farnb 

K+l 


=  0, 


hi,x=o  h(x)  /n(l-x)  y 

0  +  n+  \]  (0  +  n)2  +  4(0  +  1)(77  +  1) 


where  y  = 


2(0+1) 


(33) 


Eq.  (33)  provides  many  choices  in  the  design  of  the  manifold 
shape  with  simpler  form  than  Eq.  (26).  A  profile  such  as  h  i  (x)  = 
(Ax  +  B)~l  can  be  used  with  the  coefficients  being  determined 
by  the  substitution  to  produce: 


Fig.  5.  Profile  of  manifold  shape  for  h\(x)  =  h2(L  —  x)  and  normalized  static 
pressure  at  upper  manifold  (Z-type). 


most,  is  second  order  for  either  h\  or  h2. 


hi,x=o 

hi 


(  P  yj— - ^  P amb  ^ 

\  K*=0 


h2,x=0 

h  2 


(  r 

I  P  ,  7  P amb 

V  ,l2.x=0 


-pro  -  if 


(37) 


From  the  definition  of  M\ ,  the  static  pressure  at  the  upper  man¬ 
ifold  is: 


h\  x— o 

Pi  =  -=0/Vb 


hi 


h\,x=0 


-o-o- 


(38) 


h\(x) 


Y 

(1  -  y)x  +  y 


(34) 


The  normalized  static  pressure  at  the  upper  manifold  is: 

Px  =  fo  +  0)-  1/7(1 -O2  (35) 

h  i  hi 

The  profile  of  Eq.  (34)  together  with  the  accompanying  static 
pressure  are  shown  in  Fig.  5.  The  static  pressure  rises  and  falls 
along  the  co-ordinate  x  while  the  height  of  the  manifold  increases 
continuously. 


3.2.  U-type 


For  a  U-type  manifold,  Eq.  (23)  is  substituted  for  q  in  Eq. 
(20)  to  produce: 


Mi 

hi 


^f-pq2(\-x)2 

h2 


(36) 


Rearranging  terms  after  the  substitution  of  Eqs.  (21)  and  (22) 
for  Mi  and  M2  yields  the  second  order  manifold  equation  for 
U-type.  This  equation  appears  to  be  very  complicated  but,  at 


3.2.1.  hj(x)  =  const 

The  first  case  to  be  explored  for  the  U-type  configuration 
is  the  upper  manifold  with  a  uniform  cross-section.  The  cross- 
section  of  the  lower  manifold  is  changed  so  that  there  is  an  even 
distribution  among  the  channels.  During  the  derivation,  scaling 
variables  a ,  77  and  0  are  defined  as  in  the  previous  section  for 
Z-type  configuration  with  h 2^=1  replaced  by  /1 2,1—0,  he.: 

2 

pq 

h2,rx=o  =  a  hi,  — - -T - =  77, 

*amb™2,jt=0 

£2  h-2 

- =  6),  — —  =  h2, 

P'dmb  ^2,x=0 

[1  +  <x2n(2x  -  x2)]hj  -  (77  +  l)/z2  +  77(1  -  x)2  =  0  (39) 

The  profiles  of  the  lower  manifold  for  changes  in  the  parame¬ 
ter  a  from  0.5  to  2.0  are  given  in  Fig.  6.  For  a  =  1.0,  the  lower 
manifold  with  a  uniform  cross-section  gives  rise  to  an  even  dis¬ 
tribution  among  the  channels.  Though  this  is  the  result  obtained 
by  neglecting  the  friction  at  the  wall,  it  still  reports  that  the  U- 
type  is  more  efficient  than  Z-type  configuration  in  producing  an 
even  distribution.  Unlike  the  Z-type,  a  lower  manifold  of  the  U- 
type  is  required  to  have  an  increasing  or  decreasing  cross-section 
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Fig.  6.  Profile  of  lower  manifold  for  h  i  (x)  =  const  and  normalized  static  pressure 
at  upper  manifold  (U-type). 


in  the  direction  of  x  depending  on  its  size,  i.e.: 

?i  =  1  +  0  +  a2 T7(2x  -  x2)  (40) 

The  static  pressure  increases  continuously  along  the  co-ordinate 
x  regardless  of  the  manifold  shape.  The  Z-type  manifold  could 
have  an  increasing  or  decreasing  profile  of  static  pressure 
depending  on  its  size  and  design.  By  contrast,  the  static  pressure 
of  the  U-type  manifold  is  not  sensitive  to  the  manifold  shape  but 
is  changed  by  the  loss  or  the  gain  of  the  inertia. 


3.2.2.  Ii2(x )  =  const 

The  next  case  of  constant  h 2  also  suggests  scaling  variables  ft, 
n  and  0.  The  physical  meaning  of  each  variable  is  the  same  as 
that  discussed  above  for  the  Z-type  manifold  (see  Section  3.1.1), 
i.e.: 


,  in 

h  1  v=o  =  — , 

P 

Q 


pq 


farnb 


=  0, 


^arnb^2 
hi 


=  n. 


h\,x=0 


=  h, 


(i  +  0  +  n(2x  -  x2))h\  -  (i  +  0  +  p2n)h  1 
+  p2n{i  -  a-)2  =  0  (4i) 

This  case  produces  the  same  profile  for  the  manifold  shape  as 
that  shown  in  Fig.  6  if  a  is  replaced  by  1/p,  which  is  expected 
from  the  symmetry: 

Pi  =  jl(i  +  0  +  p2 n)  -  ^p2n(i  -  x)2 

hi 

=  1  +  0  +  TJ(2x  -  x2)  (42) 


3.2.3.  hj(x)  =  h2(x) 

The  relationship  hi(x)  =  ln{x)  =  const  is  the  only  solution 
to  satisfy  the  condition  h\ U)  =  hnix).  The  situation  where  of 
a=  1.0  in  Section  3.2.1  corresponds  to  this  case  as  shown  by 
the  data  in  Fig.  6,  i.e.: 

Pi  =  1  +  0+  I7(2x  -  x2)  (43) 

4.  Conclusions 

A  method  has  been  developed  for  designing  fuel-cell  man¬ 
ifolds  to  distribute  the  reactants  among  the  channels  within 
bipolar  plates.  Expressions  for  the  manifold  height  have  been 
derived  for  Z-type  and  U-type  configurations  after  neglecting 
the  friction  term.  Criteria  for  determining  the  validity  of  the  cur¬ 
rent  model  are  also  proposed.  Some  cases  are  explored  to  find  the 
specific  design  that  satisfies  the  derived  manifold  equation  for 
each  configuration.  The  static  pressure  within  each  manifold  is 
calculated  and  can  be  used  to  estimate  the  pressure  requirement 
at  the  inlet  to  drive  the  specific  flow  rate.  During  such  derivation, 
scaling  variables  are  defined  and  are  based  on  the  operating  con¬ 
ditions  as  well  as  on  the  geometric  configurations  of  the  channel 
and  the  manifold. 
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